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Abstract. We present simultaneous high time resolution (1-10 Hz) X-ray and optical observations 
of the persistent LMXBs ScoX-1 and V801 Ara(=4U 1636-536). In the case of ScoX-1 we find that 
the Bowen/Hell emission lags the X-ray light-curves with a light travel time of ~ 1 1 — 16s which is 
consistent with reprocessing in the donor star. We also present the detection of three correlated X- 
ray /optical bursts in V801 ara. Although this latter project is still in progress our preliminary results 
obtained by subtracting the Continuum light-curve from the Bowen/Hell data provide evidence of 
orbital phase dependent echoes from the companion star. 
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INTRODUCTION 

Low mass X-ray binaries (LMXBs) are interacting binaries containing a low mass donor 
star transferring matter onto a neutron star (NS) or black hole. Mass accretion takes place 
through an accretion disc, and with temperatures approaching ~ 10 8 K, such systems are 
strong X-ray sources. The mass transfer rate supplied by the donor star, M%, is driven by 
the binary /donor evolution and for M2 > M cr \ t ~ 10 _9 M Q yr _1 results in persistently 
bright X-ray sources. In these binaries optical emission is dominated by reprocessing 
of the powerful X-ray luminosity in the gas around the compact object which usually 
swamp the spectroscopic features of the weak companion stars. In this scenario, dynam- 
ical studies have classically been restricted to the analysis of X-ray transients during 
quiescence, where the intrinsic luminosity of the donor dominates the optical spectrum 
of the binary. 

Fortunately, this situation has recently changed thanks to the discovery of the narrow 
emission components arising from the donor star in the prototypical persistent LMXB 
Sco X-l fl 1 3M . High resolution spectroscopy revealed many narrow high-excitation emis- 
sion lines, the most prominent associated with Hell A4686 and NIII AA4634-41 / CIII 
A A 4647-50 at the core of the broad Bowen blend. In particular, the NIII lines are pow- 
ered by fluorescence resonance through cascade recombination which initially requires 
EUV seed photons of Hell Lya. These very narrow (FWHM < 50km s~ ! ) components 



move in antiphase with respect to the wings of Hell A 4686, which approximately trace 
the motion of the compact star. Both properties (narrowness and phase offset) imply 
that these components originate in the irradiated face of the donor star. We repeated this 
experiment with success in other 8 sources (see e.g. [21]) which imply that mass donor 
star irradiation and subsequent Bowen blend fluorescence is a common occurrence in 
persistent LMXBs. 



Echo-Tomography 

Echo-tomography is an indirect imaging technique which uses time delays between 
X-ray and UV/optical light-curves as a function of orbital phase in order to map the 
reprocessing sites in a binary [fl2|] . The optical light curve can be simulated by the 
convolution of the (source) X-ray light curve with a transfer function which encodes 
information about the geometry and visibility of the reprocessing regions. The transfer 
function (TF) quantifies the binary response to the irradiated flux as a function of the 
lag time and it is expected to have two main components: the accretion disc and the 
donor star. The latter is strongly dependent on the inclination angle, binary separation 
and mass ratio and, therefore, can be used to set tight constraints on these fundamental 
parameters. Successful echo-tomography experiments have been performed on several 
X-ray active LMXBs using X-ray and broad-band UV/optical light-curves. The results 
indicate that the reprocessed flux is dominated by the large contribution of the accretion 
disc (e.g. [l8ll. [fl2ll ) which dilutes the reprocessed signal arising from the companion. 
Only [9] and [4] have reported the detection of reprocessed pulsations in the X-ray 
binaries 4U 1626-67 and Her X-l respectively. More recently, [7] have reported possible 
evidence for a reprocessed signal from the companion through the detection of orbital 
phase dependent echoes in EXO 0748-676. 



BOWEN BLEND ECHO-TOMOGRAPHY OF LMXBS 

Based on our previous works where we detect reprocessed Bowen emission from the 
donor star in many LMXBs, we developed the idea of Bowen blend echo-tomography. 
To do this, we manufactured two narrow (FWHM =100 A) interference filters centered 
at the Bowen and Hell lines (A e ff=4660A) and a featureless region at A e ff=6000A. These 
filters were built to be used with ULTRACAM, a triple-beam CCD camera which uses 
two dichroics to split the light into 3 spectral ranges: Blue (< A3900), Green(AA3900- 
5400) and Red (> A5400). It uses frame transfer 1024x1024 E2V CCDs which are con- 
tinuously read out, and are capable of time resolution down to 2 milliseconds by reading 
only small selected windows (see ^\ for details). The combination of ULTRACAM and 
the narrow band filters allows us to obtain continuum- subtracted light-curves of the high 
excitation lines where the response of the donor's contribution is amplified by suppress- 
ing most of the background continuum light, which is associated with the disc. As a first 
step, we decided to undertake an echo-tomography campaign on the brightest LMXB of 
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FIGURE 1. Upper panel: Simultaneous RXTE and ULTRACAM (SDSS u' , Bowen+Hell and Contin- 
uum) observations of Sco X-1/V818 Sco obtained on the night of 18 May 2004. Lower panel: Zoom of 
W3 showing a 1000s stream of X-ray (above) and Bowen+Hell non continuum subtracted (below) data. 
The variability appears clearly correlated. 



all, Sco X-l which is also the prototypical Z — source. Here we also present preliminary 
results of a new echo-tomography campaign on the X-ray burster V801 Ara. 



Echoes from the companion in Sco X-l 

Simultaneous X-ray and optical data of Sco X-l were obtained on the nights of 
17-19 May 2004. For the X-ray observations we used the Proportional Counter Array 
(PCA) onboard the Rossi X-ray Timing Explorer (RXTE) satellite. The STANDARD- 1 
mode, with a time resolution of 0.125s, was used for the variability analysis and the 
STAND ARD-2 mode data, with a time resolution of 16s, was used for the spectral 
analysis. Sco X-l was observed during 15 RXTE windows of 16-31 minutes length 
spread over the three nights and yielding 20 ks of data. About 80% of these observations 
were also covered with simultaneous optical data. 

The optical data were obtained with ULTRACAM at the Cassegrain focus of the 4.2m 
William Herschel Telescope (WHT) at La Palma. We used our interference filters in the 
green (Bowen+Hell) and red (Continuum) arms of the intruments. A standard SDSS 
u' (A 3543 A) filter was also mounted in the blue channel. The images were reduced using 



TABLE 1. TF and ICF analysis for W3 



ICF±lc7(s) 



To±1ct(s) Ato±1ct(s) 



Bowen+Hell CS 
Bowen+Hell 

SDSS u' 
Continuum 



14.3-16.3 
11.8 ±0.1 
9.5±0.7 
9.0±0.7 



13.5±3.0 8.5±3.5 

10.75 ± 1.25 6.25 ±1.5 

9.0±1.5 5.5±2.0 

8.5 ±1.0 5.5 ±1.5 



the ULTRACAM pipeline software with bias subtraction and fiatfielding. The exposure 
time was initially set to 0. 1 s but was increased to 0.25s depending on weather conditions. 

In order to look for correlated variability we performed a comprehensive cross- 
correlation analysis by using the entire data set in blocks from 2 min to 30 min (i.e. 
one RXTE window). To do this, we used a modified version of the ICF (Interpolation 
Correlation Function) method which is explained in great detail by [6]. Significant ICF 
peaks were found only on the night of 18 May, when the system stayed at the Flaring 
Branch and the X-ray flux rose up to 1.3 x 10 4 counts s . During this night the system 
shows long episodes of flaring activity (fig. Q] upper panel) and correlated variability be- 
tween the X-ray data and the optical lightcurve was observed (fig. Q] lower panel). We 
obtain highly significant correlation peaks in the three considered filters by choosing 
different data blocks during W3 (see fig [Q), especially during the second part of this 
window where strong correlated variability is present. W3 is centered at orbital phase 
0.52, i.e. very close to the superior conjunction of the donor star, when the irradiated face 
presents its largest visibility with respect to us. In many cases we find ICF time-lags of 
about lis between the X-ray and Bowen+Hell light-curves but we finally selected for a 
detailed analysis a 3 min block of the W3 where the correlation peaks are particularly 
significant. Tabled] lists the ICF results obtained for the 3 bands considered (i.e. SDSS 
u', Bowen+Hell and Continuum). 

With the aim of amplifying the signal from the donor we subtracted the Red (Continuum) 
from the Green(Bowen+HeII) channels. However, the Continuum is also strongly corre- 
lated with the X-ray emission, and hence, the subtraction process reduces the correlation 
level of the signal and add noise. In order to determine the optimal amount of continuum 
to be subtracted, we computed several ICFs after subtracting a fraction cf (with cf in the 
range 0-1) of the red continuum from the Bowen+Hell light curve. We find that the ICF 
peak of the continuum subtracted lightcurve (Bowen+Hell CS) smoothly shifts to longer 
delays for larger cf values, as expected if more disc contribution is subtracted (see up- 
per panel in fig©. We obtain a delay in the range 14.3 — 16.3s for cf = 0.6 — 0.8, when 
~ 60% of the continuum light is subtracted. For higher cf values the ICF are nosier and 
the peaks are not significant. 

As a second step in the data analysis we fitted the ULTRACAM light-curves with 
the results of convolving the RXTE data with a set of synthetic transfer functions. For 
simplicity we have used Gaussian TFs with two free parameters: the mean delay (To) 
and the standard deviation of the delay (At): 
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FIGURE 2. Upper panel: ICF for the selected 3min block of W3. The ICFs have been computed 
using the X-ray data and a Bowen+Hell light curve obtained by subtracting from the green channel the 
Continuum data scaled by a factor cf = (i. e. no subtraction), 0.4 and 0.8 . The dotted line shows the 
3(7 confidence level and the dashed vertical lines time-lags of 11.8±0.1s, 13.3 ± 0.1s and 16.0±0.3s 
respectively. Lower panel: Expected range of delays for both the companion and the accretion disc 
considering i = 50°, M\ — 1.4M and the orbital parameters of Sco X-l. The black region shows the 
1 a delays obtained for the Bowen+Hell, and Continuum data. The shadowed region represents the range 
of delays obtained by subtracting the Continuum from the Bowen+Hell light-curve. The white dashed 
line is the most probable delay for the companion according to fioll and using a = 6°. 



where A is a normalization constant. We have applied this technique to all data sets 
and find To values completely consistent with those obtained in the cross-correlation 
analysis. A summary of the results is presented in Tabled] 

Figure [2] (lower panel) shows the range of expected time-delays for both the compan- 
ion and the accretion disc as a function of the mass ratio (q = M2/M1). The time-lags of 
~ 10 — lis obtained for the Bowen+Hell data during W3 lies at the lower edge of the 
companion region, and shows that part of this emission must arise from reprocessing on 
the donor star. Regarding the continuum subtracted data, the TF method yields a time- 
lag of 13.5 ± 3.0s which is perfectly accommodated in the companion region (fig. [21 



TABLE 2. Con-elated X-ray/optical Type I burst in V801 Ara 



BURST | TIME (MJD) q> or bital ICF-h' band(s) ICF-Bowen+HeII(s) ICF-Continuum(s) 

Burst- 1 54271.04381 0.22 1.45 ±0.02 1.73 ±0.05 1.73 ±0.05 

Burst-2 54272.09179 0.85 1.75 ±0.08 1.93 ±0.14 2.17 ±0.08 

Burst-3 54272.20536 0.57 2.51 ±0.06 2.93±0.06 2.33±0.05 



lower panel), while the ICF method gives a narrower time-lag range (14.3 — 16.3s). This 
is larger than expected for the maximum response of the companion and suggests that 
either the mass of the NS or the orbital inclination are higher than the assumed values. In 
addition the Continuum data shows delays consistent with reprocessing in the accretion 
disc whereas the u' band present larger delays. This is probably due to the presence of 
many high excitation lines arising from the companion in the spectral range covered by 
this filter. We note that correlated variability is also found during W5 (fig. CO) when the 
system shows less activity. In this case we find delays consistent with reprocessing in an 
accretion disc with a radial temperature profile, (see [1 1] for details) 



Phase-dependent echoes from the companion star in 
4U 1636-536/V801 Ara 

Encouraged by the successful echo-tomography work on Sco X-l we decided to 
undertake a new campaign on V801 Ara where reprocessed Bowen emission from the 
companion has been reported [3]. This system is also a well-known X-ray burster which 
offers the opportunity of exploiting burst events to detect reprocessed signals from the 
companion at different orbital phases. We observed V801 Ara on 19-20 June 2007 with 
ULTRACAM attached to the Very Large Telescope (VLT) at Cerro Paranal (Chile). 
We used the same instrument set up as for the Sco X-l campaign. We used a time 
resolution of 0.5s, although it was increased up to Is during the 2nd night when weather 
conditions decreased. We also obtained simultaneous RXTE data covering the ground 
based campaing. Here we present a preliminary ULTRACAM-pipeline data reduction 
together with the STANDARD- 1 RXTE data set at T exp = 0. 125s. 
Three simultaneous X-ray/optical Burst were detected in the course of the 2 nights (see 
fig E] left panel). Table [2] lists the times of the bursts and the orbital phase at which the 
events occured according to the ephemerides in |3|]. 

ICF functions were computed in the same way as in Sco X- 1 by selecting 200s blocks 
of data around each burst. We find in all the cases very clear ICF peaks centered at 
positives delays (i.e optical event occurs after X-ray event). Figure[3]shows the X-ray and 
Bowen+Hell light-curves during Burst- 1 together with the resulting ICF function (right 
panel) whereas table [2] lists the obtained time lags for all bands considered. We measure 
longer delays during Burst-3, at <p or & ~ 0.57, when the visibilily of the companion is 
larger. Furthermore, similar as for Sco X-l, the delay obtained for the Bowen+Hell 
filter is longer than those observed in the other bands. 

In order to amplify the signal from the donor we subtracted the Continuum emission 
scaled by a continuum factor (cf) from the Bowen+Hell lightcurve. Figure 0] shows the 
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FIGURE 3. Left panel:40s streams of X-ray (above) and Bowen+Hell (below) data around Burst- 1 . The 
smearing in the shape of the optical event is clearly seen in the figure. Righ panel: ICF obtained through 
the left panel lightcurves. A clear peak centered at 1 .73 ± 0.05 is found. 



delays obtained for each burst vs cf. We also mark (grey regions) the expected range 
of delays for the companion (following iflQIl ) by considering an orientative set of orbital 
paramenters and the orbital phase at which the events were detected. In the three cases 
the time-lags measured are consistent with reprocessing on the companion star. 



CONCLUSIONS 

Following the discovery of narrow high excitation emission lines arising from the com- 
panion star in many LMXBs we are undertaking several Bowen blend echo-tomography 
experiments using ULTRACAM+narrow band filters and RXTE. Our first work on the 
prototypical LMXB Sco X-l presents the detection of echoes consistent with reprocess- 
ing on the companion star. Correlated variability was found when the system stayed at 
the Flaring Branch and displayed high and low frequency variability [11]. Preliminary 
results of a new campaign on the X-ray burster V801 Ara yield evidence of orbital- 
phase dependent echoes from the companion during the three burst which were detected. 
This latter work is still in progress and we are currently refining our results by fitting 
time-delayed transfer functions to our data. More observations of other targets with this 
promising technique are currently underway. 
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FIGURE 4. Time-lag vs Continuum subtraction factor for the three burst detected. The grey regions 
show the expected range of delays for the companion assuming q = 0.265 ([3]), M^s = 1.4M Q , disc 
flaring angle (a) = 6°, orbital inclination in the range 40-60° and the orbital phase corresponding to each 
burst. 
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